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Abstract: Wo used a complete tight-binding band structure of graphene 
nanoribbon to obtain, for the first time, analytical techniques for observ- 
ing photon assisted transport, and dynamic localization of electrons in the 
graphene nanoribbons. When the ribbons are subject to a multi-frequency 
dc-ac field, photon assisted replicas show up at rather strong drive force. 
The strong dependence of the photon peaks on ac amplitudes allow for 
high-harmonic dynamic oscillations at these amplitudes. We identified re- 
gions of positive differential conductivity where a nanoelectronic graphene 
device may be operated as a small signal amplifier. Our research has also re- 
veal another quantum mechanical phenomenon, fractional photon assisted 
transport, when the stark factor r > 1. 
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1. Introduction 

Photon assisted processes in systems are quantum phenomenon in the pres- 
ence of high electric fields in which electrons can absorb or emit one or several 
photons. These absorption or emission of photons are seen as new conduction 
channels in the system. The observation of photon assisted process in solid sate 
materials is central ingredient for terahertz (THz) generations. In semiconductor 
nanostructures, photon assisted peaks (PAPs) has emerged as a powerful tool 
for investigating quantum transport phenomenon. PAPs up to 7 THz photons 
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has been reported in [1] and can persist up to room temperature. Recent studies 
have addressed electronic properties of confined graphene structures hke dots, 
nanowires, nanotubes (NTs) and nanoribbons (NRs). In particular, nanoribbons 
have been suggested as potential candidate for replacing electronic components 
in future nanoelectronics devices [2, 3]. Research have also shown that graphene 
is a suitable candidate to examine photon assisted processes in Dirac systems 
[4]. For this reason among others, photon assisted tunneling was studied quite 
recently in graphene [5] and also [6]. In NTs, PAPs occurs because the current 
dynamics in the system can be seen as Bloch oscillations when an ac field is 
added to dc field but with frequency somewhat below the Bloch frequency as- 
sociated with the dc field. The large ac amplitude opens up transport channels 
which can be seen as photon peaks. Regions of positive differential conductiv- 
ity at certain ranges of the dc field are formed. As a consequence, domainless 
amplification of THz frequencies can be achieve [7] . 

When a superlattice is irradiated with an intense dc-ac field (laser pulse), it can 
make charged particles localized around the initially incident region, this results 
in what is called dynamic localization [9, 10]. The combined effect of the dc-ac 
field is crucial and can cause destructive and constructive effects in electron 
dynamics. In the latter, regions of positive dynamic conductivity are formed. 
The destructive case is manifested by what the dynamic localization unless the 
ratio of Bloch frequency to the applied frequency is not an integer and that the 
argument of a Bessel function whose order is an integer is not a root [8] . 

In this paper we investigate the possibility of photon assisted processes (replicas) 
in GNRs using Boltzmann transport equation in the relaxation time approxima- 
tion. Motivated by the fact that photon assisted tunneling has recently appeared 
in [5] where a very simple Dirac spectrum describing only low energy electrons 
was used, we explore not the tunneling of Dirac electrons through barriers but 
the operation regions for possible THz amplification using full tight binding 
spectrum of graphene. We will show that the I-V characteristics demonstrate 
regions of positive differential conductivity (PDC) where a graphene nanodevice 
can be operated as a small signal amplifier. We noted that a dynamical response 
of the electrons to terahertz radiation can depend on details of scattering pro- 
cesses, this effect could lead to the relaxation towards a thermal equilibrium 
which we relax here. Again we ignore transient processes, that is we assume 
that t»T [13]. 

The rest of the article is organized as follows; we derived the photon assisted 
current density in section 2. Using the photon currents, we study the photon 
replicas and dynamic localization in the presence of monochromatic and bi- 
harmonic frequencies in section 3. I-V characteristic graphs are plotted and 
discussed in details in section 4 and the paper finally concludes in section 5 
outlying some further investigations. 
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2. Photon assisted current density 

The current density of electrons in GNR subject to a general multi-frequency 
field, E{t) = Eo + EjCos{ujjt + aj) is [12] 
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where J„(a;) are Bessel functions of n*^ order. /3j — elEj/hui, /3o = elEo/H and 
Jo is the peak current define as 
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M = Ti/iM + 1) and I = ^/3a/2 for aGNR, M ^ [2 + l/s)Ti/{J\f + 1) and 
I = a/2 for zGNR. In order to have more simplified equation that describes the 
phenomenon we are considering, we adopt the following definitions; 
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with ujjT = E* /Ecr- E* is the electric field at which an electron emits or absorbs 
a photon when displaced by quasi-lattice period I. Ecr is however, the critical 
field at which carriers get their peak velocity. With these choice of variables 
we make the replacement /3or + rijUjjT — > E^/Ecr + rijE* / Ecr, so that Eq.(2.1) 
becomes 
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where f{x) are superposition of weighted photon replicas of pure dc differential 
current density. If the above equation is compared to the pure dc current density 
in reference [12], we see that the photon peaks are displaced by multiples of 
E* / Ecr with amplitudes J^. Note that if r = 1 and Uj = 1, the Bloch frequency, 
ujb {= elEo/h) associated with the static field will coincide with the applied 
frequency loj and no photon replicas can be seen. Therefore it is a requirement 
that r 7^ Tij and Uj 7^ 1, if PAPs are to be observed. In the following section, we 
will consider two cases of Eq.(2.2); the monochromatic and bichromatic cases 
in the presence of the static electric field. 
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3. Large-signal dynamic localization and photon-assisted 
replicas 

3.1. Graphene nanoribbon in monochromatic field 

When an ac field of frequency, lo is applied to the GNR, we have n — 1, ujjyi ~ 
and oij = in Eq.(2.2) whieh becomes 
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Ecoslvut). 
(3.1) 

The coefficient of Ecos{ujt) in Eq.(3.2) is an ac part of the drive field. It is seen 
as a large-signal dynamic conductivity at the derive harmonic frequency lo. The 
fundamental derive frequency is obtained by setting u ~ ±1. To simplify the 
preceding equation further, we take v = Q to get cos{i'ujt) = 1, so that 



Jn[r-fT-]f[ r—+n— 



n— — oo 



E* 



E, 



cr 



E, 



cr 



(3.2) 



3.2. Graphene nanoribbon in bichromatic field 



Now, for an applied field of frequencies, uji and u!2, we have n = 2 and <jjj>2 = 0- 
These set ajyi = 0. One can put vi ^ and 7^ allowing Eq.(3.1) look 
like 
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We arrived at the last equation by allowing the two frequencies to be commensu- 
rable, i.e li^ilwi = |z^2|w2, and periodic, i.e with W2 = /uwi. You can see reference 
[13] and references there in for the ease of non-eommensurable frequencies. Fur- 
ther, if both 1^1 = and V2 = ^ then a more simplified equation 
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is obtained. A significant difference between the above two equations is the 
phase factor cos{v2a). Of course, among other differences, Eq.(3.4) will be an 
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inversion of Eq.(3.3) for a £ [tt/2,3tt/2]. The advantage of the biharmonic field 
over monoharmonic field is that in the former, new local structures in the I- 
V characteristics become most pronounce. We shall see this in the following 
section. 

4. Results and Discussions 

The graphs in Fig.l are j — Eq plot of Eq.(3.1) which demonstrate appearance 
of positive differential conductivity (PDC) regions at high applied frequencies, 
OJT ^ 10. Particularly in AGNR, when the Bloch frequency associated with 
the bias field falls within the ranges 9.15r"^ < ujb < 10.9t~\ 19.25t~^ < 
ub < 20.65r~^ etc., NDC is greatly suppressed and the device can be operated 
effectively. This means the system is free from space charge instability. Any small 
signal passing through the ribbon no more suffer from electrical domains and can 
thus easily get modified. The scenario is the same for ZGNR except that some 
operation points are hidden at large Eq/E* values. From the figure, because the 




Fig 1. I-V characteristic showing regions of PDC for (left) armchair and (right) zigzag rib- 
bons. ujT = 10, E = E* . 

peaks show up at Ecr/E* + n > Eq/E* and the minima at Eq/E* > n, where 
n is the number of peaks, one can deduce a general range of bias field for best 
device operation to lie within n < Eq/E* < ti + Ecr/E* or 

nwT < ujbt < nujT + 1. 

The behavior of Eq.(3.2) is shown in Fig. 2. PDC effect is observed at high even 
and odd harmonics. We chose ly = 2, v = 3 and cos{i'a) = 1. The two graphs 
are distinct, they are mirror reflections of each other. NDC in AGNR becomes 
PDC in ZGNR and vice-versa. It is clear from the two curves that depending 
on the sign of cos(^a), the operation ranges in one ribbon can be greater than 
the other. 
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To see the existence of dynamic localization in graphene nanoribbons, we plot 
Eq.(3.2) for normalized conductivity versus E / E* for default parameters lot — 5 
and ujbt = 4 in Fig. 3. 
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Fig 3. Dynamic localization states induced by ac-dc fields occur at E/E* = 3.8, 7.02, 10.09. 
u = 



Photon assisted peaks are shown in Fig. 4 (blue) for aGNR. Multi-photon res- 
onances appear when vlubt = nujT, for an integer n. If n/r is the order of the 
Besscl functions with root coinciding with one of the E/E* values in Fig. 3, a 
dynamic localization is seen. Depending on how the ratio E/E* is chosen, the 
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localized charged carriers can extend the localization beyond its vicinity and 
thus affecting neighboring centers. This might account for decrease in neigh- 
boring peak heights. The suppression at the first photon assisted peak of the 
red curve is due to the combined destructive effect of ac-dc fields at Eq = E* . 



01 234 "01 234 

Fig 4. (blue) Photon assisted peaks for E/ E* = 2.0. The peaks occur at uji, = uj, 2lj, 3uj, . . . 
(red) Dynamic localization at first peak. We chose E/E* = 7.0. 




To observe high-harmonic dynamic localization in graphene, we plot the high 
dynamic conductivity equation with E/E* for = 2 in Fig. 5. We obtained 
the localized states at E/E* ~ 5.05, 8.45, 11.65. Localized sates appear later in 
high-harmonic case compared to the normal dynamic localization. This might be 
an an indication that at high harmonics dynamic localization could disappear 
giving way for proper direct signal rectification. Fig. 6 shows photon assisted 
resonances and high-harmonic dynamic localization at the first photon assisted 
peak. 

A very important phenomenon of fractional photon assisted processes is ob- 
tained in graphene nanoribbon when the stark component r > 1, even though 
the integer photon assisted peaks are still more pronounced. In Fig. 7 fractional 
peaks arc formed at n/r = 0.5, 1.0, 1.5, 2.0, . . .. 

5. Conclusion 

We have used the complete tight-binding spectrum of graphene to show the 
phenomenon of photon assisted processes and dynamic localization in graphene. 
A signature which may be responsible for rectification and amplification of ac 
fields. We found yet another quantum mechanical behavior of graphene nanorib- 
bons, fractional photon assisted process. We have suggested the use of our theo- 
retical approach for studying terahertz generation and small signal amplification 
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Fig 5. High dynamic localization states induced by ac-dc fields occur at E/E* = 
5.05,8.45,11.65. v = 2 




Fig 6. (blue) Photon assisted peaks for E/E* = 2.5. The peaks occur at ljb = i^, 2a;, 3uj, . . . 
(red) High- harmonic dynamic localization at first peak. We chose E/E* = 5.0. 



in graphene nanoelectronic devices. 
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Fig 7. Fractional photon assisted process in graphene nanoribbon. 



References 

[I] S. Zeuner , S. J. Allen, K. D. Maranowski and A. C. Gossard, Appl. Phys. 
Lett. 69 2689 (1996). 

[2] A. K. Geim and K. S. Novoselov, Nature Materials, v.6, 183 (2007). 
[3] J. W. Son, M. L. Cohen and S. G. Louie, Nature, v.444, 347 (2006). 
[4] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov andA. K. 

Geim, Reviews of Modern Physics, v.81, 109 (2009). 
[5] A. lurov, G. Gumbs, O. Roslyak and D. Huang. J. Phys.: Condons. Matter 

24 015303, (2011), arxiv:1107.5779vl (2011). 
[6] B. Trauzettel B., Ya. M. Blanter, and A. F. Morpurgo, arXiv:0606505v2 

(2006). 

[7] S. S. Seidu, Negative Differential Conductivity and Terahertz Gnereation in 
Carbon Nanotubes, PhD Thesis, University of Cape coast. Department of 
Physics, Ghana, pp 146 (2011). 

[8] W. X. Yan, S. Q. Bao, and X. G. Zhao, arXiv:cond-mat/9801012vl (19198). 

[9] A. W. Ghosh, A. V. Kutznetsov, and J. W. Wilkins, Phys. Rev. Lett. 79, 
3494 (1997). 

[10] X. G. Zhao, Phys. Lett. A 155, 299 (1991). 

[II] H. Kroemer, arXiv:0009311 (2000). 

[12] M. Musah, S. Y. Mensah and S. S. Abukari, arXiv:submit/0418298 [cond- 

mat.mtrl-sci] (2012) 
[13] Yu. A. Ramanov, J. Yu. Romanova, and G. Morouk, arXiv:0901.0961 

(2008). 

[14] T. Hyart , A. Natalia, V. Alexeeva, J. Mattas, K. N. Alekseev, Elsivier, 
Microelectronics Journal 40 719721 (2009). 



Rabiu Musah et al. /Photon Assisted Peaks in graphene nanoribbons 10 

[15] T. Hyart, K. N. Alekseev, E. V. Thuneberg, Phys. Rev. B 77 165330 (2008). 



